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ABSTRACT 

This  report  studies  the  existence  of  non-constant  solutions  of  certain 
two-point  boundary  value  problems  for  3-component  systems  with  a  small 

,  t 

parameter  e,  under  homogeneous  Neumann  conditions  at  the  boundaries.  This 
problem  is  related  to  the  analysis  of  segregation  patterns  in  population 
models  of  3-competing  and  spatially  dispersing  species.  It  is  shown  that  the 

*"  l  ' 

reduced  problem  *  0)  has  many  non-constant  solutions  exhibiting  spatial 
segregation.  Only  a  few  of  these ,  however,  can  serve  as  valid  lowest-order 
approximations  to  solutions  of  the  original  problem  when  4  is  non-zero  but 
small.  A  singular  perturbation  construe '-ion  clarifies  which  are  in  this 
category.  The  results  of  numerical  computations  of  solutions  are  also 
illustrated. 
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A  3-COMPONENT  SYSTEM  OF  COMPETITION  AND  DIFFUSION 


Masayasu  Mimura  and  Paul  C.  Fife 


1.  Introduction 

He  consider  populations  of  N  species  S^,  S2,...,SN  in  a  bounded  habitat,  and  assume 
that  the  distribution  of  the  populations  are  determined  by  competition  of  Lotka-Volterra- 
Gause  type  and  simple  diffusion.  Suppose  U|(t,x)  (i= 1 , 2, . . . ,N )  is  the  population  density 
of  the  species  S^  ( i=1 , 2, . . . ,N) .  The  spatial  domain  is  taken  to  be  the  one-dimensional 
interval  (0,1).  *nien  we  have  the  following  reaction-diffusion  equations  governing  the 
evolution  of  the  u^: 

3u.  3  u  N 

M)  =  di  ~~T  +  (ri  "  ^  “iiVV  >  °*  *  e  (i=1.2 . N) 

3x  j=1  -*  J 


where  d^,  r^  and  a^  ( i, j=1 ,2, . . . ,N )  are  non-negative  constants.  In  ecological 
terms,  ri  is  the  intrinsic  growth  rate  of  S^  a£i  is  a  measure  of  intraspecific 
competition  of  S^,  and  a^  (i^j)  is  a  measure  of  interspecific  competition  between  the 
species.  The  boundary  and  initial  conditions  are  taken  to  be 


(2) 

and 

(3) 


3x 


( t , x )  =0  for  t  >  0,  x  =  0,1,  (i-1,2. 


,N) 


(0,x)  =  ui(J(x)  >  0  for  x  e  [0,1] 


The  case  of  two  competing  species  (N=2) 

3u 


(4) 


3t 

3u 


i  ®2ui 

■ =  di  zr +  <ri 


Accession  For 

I  FT IS  GRAAI 
j  r/i'IC  TAR 
|  r-mnnonno  •  ' 

!  J-,  If;  _ 

!  pv _ _ _ 

t  Di  si  r:  ••••.■! i-n/ 
Av.-i : ' il  liv  C 


_ 2 

V*t  ’  "2  3x2  +  <r2 


a  u  —  a  u  )u 
11  1  12  21 


a21U1  -  a22U2,U2 


11 


r  i-'il 


A 


has  been  studied  extensively.  When  r^  and  a^  (i,j=1,2)  are  fixed  arbitrarily  except 
that  the  inequalities 
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are  not  both  satisfied, 
solution  (for  instance, 
large  tine.  Furthermore 


!n  <  !i  <  In 
a21  r2  *22 

the  solution  of  (4),  (2),  (3)  tends  to  a  constant  steady  state 
[2]).  That  is,  there  never  appears  any  spatial  segregation  for 
,  in  the  remaining  case 

*11  <  r1  f  *12 
*21  r2  *22 


Kiahlmoto  [5]  showed  that  even  if  there  exist  non-constant  nonnegative  steady  state 
solutions,  these  are  unstable.  This  result  indicates  that  the  simple  dlf fusion-competition 
model  (4)  cannot  explain  spatially  segregated  distributions  of  two  competing  species,  nie 
situation  is  different  when  one  introduces  the  effect  of  cross-population  pressure  [13] 
into  (4).  It  is  shown  in  [11]  that  there  occur  non-constant  steady  state  solutions 
exhibiting  spatial  segregation  on  the  basis  of  cross-dif fusion-induced  instability  (see 
also  [10]  ). 

These  results  motivate  us  to  study  whether  or  not  systems  with  more  than  two 
components  (N  >  3)  can  exhibit  spatial  segregation.  In  this  connection,  Evans  [3]  and 
Kishimoto  [6]  have  already  presented  some  examples  of  competition  and/or  prey-predator 
systems  in  the  case  N  -  3  which  show  diffusion-induced  instability.  This  phenomenon  has 
usually  been  associated  with  activator-inhibitor  systems,  but  competitive  interaction  does 
not  fall  into  that  category.  A  few  examples  of  systems  with  stable  non-constant 
bifurcating  steady  state  solutions  were  given  in  [6]. 

Furthermore,  Kishimoto,  Mimura  and  Yoshida  [7]  have  recently  studied  the  system  (1) 
with  N  »  4  and  have  shown  that  there  appear  time-periodic,  spatially  non-constant 
solutions. 

In  this  paper,  we  will  be  interested  in  the  stationary  problem  (1)  and  (2)  with  N  ■  3 
under  the  conditions  that  two  of  the  diffusion  coefficients  are  sufficiently  small.  We 
show  the  existence  of  non-constant,  nonnegative  solutions  by  a  singular  perturbation 
construction. 
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2.  Ecological  background  and  assumption 

Uncovering  the  mechanisms  for  spatial  patterning  of  ecological  communities  has  been  a 
very  active  line  or  research  in  population  hiology.  Most  often,  spatial  diversity  in 
connected  habitats  has  been  assumed  to  be  linked  to  some  heterogeneity  in  the  environment, 
but  it  has  also  been  recognized  (Levin  [9] ,  for  example)  that  stable  patterns  can  exist  in 
otherwise  homogeneous  environments,  A  celebrated  mathematical  model  displaying  such 
features  in  the  analogous  context  of  chemical  interaction  (and  diffusion)  was  presented  by 
Turing  (15]  in  1954,  and  model  systems  of  this  general  type  have  been  elaborated 
extensively  since  then. 

Our  context,  of  course,  is  a  continuous  space-time  population  model  incorporating 
species  interaction  and  dispersal.  Interactions  of  prey-predator  type  are  analogous  to 
those  of  "activator-inhibitor"  type  occurring  in  other  disciplines,  and  are  most 
reminiscent  of  the  type  of  chemical  interaction  that  Turing  postulated.  And,  in  fact, 
small  aii^litude  patterns  for  two-species  prey-predator  situations  with  spatial  dispersal 
modelled  by  diffusion  operators  can,  in  many  cases,  be  constructed  by  standard  bifurcation 
techniques. 

Con^etition-type  interactions  are  not  of  activator-inhibitor  type,  however,  and  so  it 
is  not  too  surprising  that  bifurcation  techniques  provide  no  stable  patterns  for  two 
competing  species  models.  More  than  this  is  true.  Recently  Kishimoto  and  Weinberger  [8] 
have  shown  that  no  non-constant  stable  solutions  of  the  corresponding  boundary  value 
problem  (bounded  and  convex  domain  in  one  or  more  dimensions)  with  no-flux  boundary 
condition  can  exist.  This  raises  the  question  as  to  whether  patterns  are  possible  in 
systems  with  more  than  two  competing  species. 

The  present  paper  shows  that  patterns  are  possible  in  three-species  competition 
systems  with  diffusion.  The  technique  used  is  not  bifurcation,  but  rather  formal  singular 
perturbation  layer  analysis,  utilizing  large  differences  in  the  diffusion  rates. 

In  practical  terms,  these  results  simply  serve  to  point  out  yet  another  mechanism  for 
stable  spatial  diversity  in  homogeneous  environments.  There  has  been,  of  course,  no 
experimental  or  observational  studies  which  clearly  demonstrate  the  action  of  this 


mechanism  in  nature,  although  segregation  phenomena  in  communities  of  three  species, 
possible  competing,  have  been  studied.  For  example,  TrSgardh  [14]  studies  the  interaction 
of  three  (and  more)  species  of  pine  beetles.  When  Blastophagus  pinlperda  and  B.  minor  were 
present  on  a  tree,  their  ranges  were  observed  to  be  segregated:  B.  plniperda  occupied  the 
lower,  and  minor  the  upper  region  of  the  trunk.  But  on  a  small  island  B.  minor  was  not 
seen,  and  the  range  of  the  other  species  was  enlarged.  A  third  species,  Ips  lonqicollls, 
living  more  in  the  outer  part  of  the  bark,  coexisted  throughout  the  ranges  of  both 
Blastophagus  species.  Crombie  [1]  observed  a  similar  ecosystem  with  Khizoperta  dominlca, 
Sltotroga  cerealella,  and  Orvzaephilus  surinamensis.  Competition  among  three  or  more 
species,  irrespective  of  spatial  distribution  considerations,  is  recently  under 
investigation  by  many  researchers:  for  a  survey  of  some  of  this  work,  see  Schoener  [12). 

We  shall  investigate  three-species  competition  systems,  and  shall  look  for  spatial 
distributions  similar  to  those  mentioned  above,  observed  by  Trlgardh  and  Crombie. 


Consider  the  stationary  version  of  (1)  and  (2)  with  N  »  3, 


d2u. 


0  =  d 


1  T  +  (rl  ‘  *11U1  ‘  *12U2  “  al3U3)U1 
dx 


(5) 


\  0  "  <*2  ~T  +  <r2  -  a21U1  •  *22U2  *  a23u3)u2  x  e  (°'1>  ' 

3x 


d2u. 


°  =  d3  ~T  +  ,r3  ‘  a31U1  '  “32U2  "  a33U3>U3 
V  dx 


and  the  boundary  conditions 
(6) 


du. 


dx 


0,  x  »  0,1,  (i-1,2,3) 


(Hi) 


and 

aij  (i,} 

23  <  _2 

*32 

a33  r3 

In  ecological  terms,  this  implies  that  when  is  absent,  Sj  can  survive  and  Sj 

becomes  extinct  in  competition.  That  is,  when  u^  =  0, 


lim  (u.(t,x),  u.(t,x))  “  (“ — 
t-H»  1  a22 


,0) 


( de  Mottoni  [2  ]  ) . 


(H2) 


-  —  <  -  and  -  <  —  <  -  , 

a22  <  r2  a21  a33  r3  a31 

which  imply  that  in  the  absence  of  S3  (resp.  S2>,  S1  and  S2  (resp.  Sj)  may  coexist  in 
competition.  (See  Figure  1). 


Remark  1 .  The  argument  In  this  paper  will  also  handle  the  case  when  S2  survives  to  the 
exclusion  of  S^,  when  S3  is  absent/  the  other  assumptions  above  remaining  unchanged. 

^H3j  and  S3  are  both  slowly  diffusing  , 

2  2 

that  is,  d1  =  1,  d2  =  e  and  d3  =  de 
for  some  constant  d  >  0  (0  <  e  <<  1)  . 

We  will  impose  some  conditions  on  r^  and  a^  (i,j=t,2,3)  in  addition  to  (HI)  and 
(H2)  later.  Note  that  the  question  as  to  whether  all  three  species  can  coexist  in  the 
absence  of  diffusion  is  not  relevant  to  the  present  study. 
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where  £  <u.|(12\3)  <  £  <  u^(13\2))  is  arbitrarily  fixed,  and  consider  the  probiem 

A, 

dx 

subject  to  the  boundary  condition  (8). 

Theorem  1 .  Fix  £  arbitrarily  so  that  maxtu^ ( 13\ 2 ) ,  j  u(12\3>)  <  £  <  u^(12\3).  Then  the 
problem  (10),  (8)  has  a  non-constant  nonnegative  solution  v^lx/C)  satisfying 

u  1  ( 1 3\ 2 )  <  v1(x;£)  <^(12X3)  . 


Proof ;  See  Mimura  [10], 

Using  this  function  v^(xj£),  we  obtain  a  solution  of  the  problem  (7), 
follows: 


(ID 


v*(x;£)  , 

r2  "  S21 V1 ^  * 
a22 

0 


r3  "  a31V1(X,5) 


*33 


(v1  >  £ ) 

(v1  <  £)  , 

(v1  >  £ ) 

(v1  <  £)  (Figure  3)  . 


(8)  as 


The  solution  (11)  shows  that  S2  and  S3  are  coexisting  with  spatial  structure  of 
segregation.  This  is  an  interesting  phenomenon  because,  in  the  absence  of  Sj,  S3  always 
becomes  extinct  under  the  assumption  (HI).  That  is,  the  presence  of  Sj  makes  the 
coexistence  of  S2  and  S3  possible. 


v2(x;£)  v3(x;£) 

Figure  3A 


v,(x;£) 


Figure  3B.  "Outer"  solution 


We  next  consider  the  problem  of  whether  the  solution  of  the  reduced  problem  (11) 
becomes  a  lowest  order  approximation  to  a  solution  of  the  original  problem  (5)  when  e  is 
sufficiently  small  but  non-zero.  The  solution  will  be  sought  in  separate  regions  named 
"inner"  and  "outer"  regions  and  then  matched  appropriately.  It  is  expected  that  the 
function  (11)  becomes  an  "outer"  solution  outside  a  neighborhood  of  x  =  x*,  where  x*  i 
a  point  where  v^x*)  =  £.  Therefore  we  seek  here  an  "inner"  solution  in  the  neighbor¬ 
hood  of  x  =  x*.  We  use  the  transformation  z  =  (x-x*)/e  in  (5)  to  obtain 


(12) 


d2u 


0  ’  "T1  +  e2(r1  "  a11U1  "  *12U2  •  a13U3>U1  ' 
dz 

d2u2 

0  -  —  *  (r2  -  *21ui  “  a22U2  '  a23U3>U2  ' 
dz 

d2u3 

V  °  ”  d  7T  +  <r3  ‘  *31U1  '  a32U2  -  a33U3)U3  * 

v  dz 


By  putting  e  -  0  and  then  u^tz)  •  v^x*)  -  ,  we  obtain  a  2-component  system  with 
respect  to  u-  end  u3,  which  is  an  approximation  to  (12): 

d2u 

o  -  — +  (R2(5)  -  a22U2  -  a2303>u2 
dz 


(13) 


A, 

0  -  d  — t  (R3(C)  -  a32u2  -  a33u3)03 
V  dz 


z  e  r  , 


where  R^(£)  *  r  -  a^C  (i  “  2,3).  The  boundary  conditions  at  z  “  ±*>  are  specified  as 

r 


(14) 


R,(C) 
0,(  — )  -  -r — 
a22 


03(— )  -  0 


o2(+->  “  0 


and  \ 


R,(5> 

„  (+.>  .  -A -  . 

v  33 


tm 


We  must  study  the  existence  of  a  nonnegative  solution  (U2(z),  U3(z)>  to  the  problem  (13), 
(14). 

Here  we  assume 

(H5) 


*21  >  T_2 
*31  r3 


which  implies  that  the  surfaces  r2  -  a21Uj  -  a22u2  -  a23u3  “  0  and  r3  ”  a31u1  ”  a32u2  ” 
a33u3  *  0  intersect  in  R2  »  { (u^ ,u2 ,u3 ) |uA  >  0  (i  »  1,2,3)}.  That  is,  two  lines 

Rj(S>  “  a22u2  "  *23U3  "  0  and  R3*^)  “  a32U2  "  a33U3  "  °  intersect  in 

*+  ”  ^u2'u3^u2  *  ®*  u3  >  tor  each  5  satisfying  ^(23x2)  <  £  <  ^(23X3).  Thus,  we 
find  that  as  ?  increases,  the  kinetic  system  of  (13)  takes  the  form  of  a  monostable 
system  (0  <  5  <  u^(23\2)),  a  bistable  one  (u^(23X2)  <  £  <  ^(23X3))  and  then  a 


monostable  one  (^(23X3)  <  C)  again. 
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Theorem  2.  There  exists  a  value  5*  satisfying  u1<23  2)  <  £*  <  Uj(23  3)  such  that  the 
problem  (13),  (14)  has  a  solution  (O^fzjC*),  1)3(1:  £*))•  (Figure  4). 


u 


Figure  4. "Inner"  solution 

Proof;  This  is  a  consequence  of  Theorem  a2  in  the  Appendix,  with  the  following  notational 
changes:  subscripts  1  in  that  theorem  are  to  be  replaced  by  3  for  the  present 
context:  r  1  is  to  be  Rj/d,  r2  is  Rj,  and  a^  is  *3 j/d* 

Remark  1.  Theorem  2  does  not  show  explicitly  the  value  of  5*.  However,  when  a22  “  “33 
and  a32  “  a23,  then  5*  *  (r2  -  r3>/<a2i  “  a31 >  (Gardner  [4]). 

Remark  2.  The  uniqueness  of  such  a  5*  is  not  known,  but  has  been  confirmed  numerically. 
Thus,  when  we  assume 

(H6)  (^(23X2),  u1  (23\3 ) )  C  (^(13X2),  ^(12X3))  , 

we  obtain  "outer"  solutions  (v^xjC*),  v2(x;!i*),  Vj(x;£*)),  as  well  as  "inner"  solutions 

x—  X*  X”  X* 

(Dj (x*» £*) ,  U2(— - — ;5*),  U3 (— - — :£*)),  where  5*  is  determined  by  Theorem  2  and  x*  is 
found  from  v. (x*:5*)  “  5*. 


The  final  task  is  to  obtain  a  solution  of  the  original  problem  (5),  (6)  when  c  is 
sufficiently  small  but  nonzero  which  corresponds  to  the  "outer"  solution  of  (7),  (8)  and 
the  "inner"  solution  of  (13),  (14)  constructed  above.  Unfortunately  we  have  not  yet  been 
able  to  prove  the  existence  of  this  solution.  We  show,  however,  numerical  experiments  to 
supplement  the  analysis  given  in  this  paper,  thus  further  confirming  the  existence  of 
solutions  exhibiting  spatial  segregation  between  two  of  the  species  (Figure  5). 


Figure  5 


4.  Concluding  remark 


We  have  discussed  non-constant  steady  state  solutions  of  a  three  component  system  of 
competition  and  diffusion.  In  this  system,  we  have  assumed  that  species  1  and  2,  as  well 
as  1  and  3,  can  coexist,  but  2  and  3  are  of  exclusive  interaction,  that  is,  either  2  or  3 
becomes  extinct.  <nius,  a  singular  perturbation  technique  constructs  non-constant  steady 
state  solutions  exhibiting  a  striking  segregated  pattern  in  the  species  when  species  1 
diffuses  fast  but  2  and  3  diffuse  slowly.  The  argument  presented  here  is  also  valid  in  the 
situation  where  species  2  and  3  coexist,  but  2  eliminates  1  as  well  as  3. 


Appendix 


Consider  the  ays tee 


(A1) 

where 


(u  +  uf(u,v)  »  0 
“ 

Vxx  *  vg(u'v)  “  0  ' 


f  -  r,  -  anu  -  a12v  end  g  -  r2  -  a21u  -  a22v  . 

(Ail  of  the  following  extends  easily  to  a  much  wider  class  of  functions  (f,g).  He  shall 
not  pursue  that  generality.) 

Definition.  The  pair  (f,g)  of  this  fore  with  r^  >  0  and  a^  >  0  (i,j  «  1,2)  is  said 
to  belong  to  class  a  if  r2/r^  >  *22^*12* 

For  (f,g)  e  a,  we  shall  be  concerned  with  conditions  under  which  there  exists  a 
solution  of  ( A1 )  for  x  e  It  satisfying 

<A2>  (u,v)(-«)  ■  (0,r2/a22),  (u,v)(+«)  »  (r^/a^.O) 

With  no  loss  of  generality  we  shall  assume,  for  simplicity,  that  r2  ”  a22  and  r1  =  a^. 
This  may  always  be  achieved  by  rescaling  u  and  v.  Written  as  a  system  of  4  first  order 
equations,  (Al)  becomes 

(  • 
u  "  p 

Ip-  -uf(u,v) 

(A3)  \  • 

V-q 

^  q  -  -vg(u,v) 

with  (f,g)  e  a,  where  •  ■  .  Solutions  of  the  system  will  be  denoted  by  X(x)  where 

X  ”  (u,p,v,q) .  Clearly  the  points  A  -  (0,0, 1,0)  and  B  -  (1,0, 0,0)  are  rest  points. 

Our  main  question  will  be  whether  there  exist  trajectories  connecting  them.  The  projection 
of  X(x)  onto  the  (u,v)  plane  will  be  denoted  by  U(x),  where  U  »  (u,v). 

Proposition  1.  The  rest  point  A  has  a  local  two-dimensional  unstable  manifold  M  which 
may  be  parameterized  by  points  in  the  (u,v)  plane  in  a  neighborhood  of  (0,1).  More 
precisely,  to  each  such  point  (u,v)  f  (0,1),  there  corresponds  exactly  one  trajectory 
on  M,  and  there  exist  no  other  trajectories  issuing  from  A.  The  trajectories  themselves 
may  be  parameterized  by  a  single  real  parameter. 

Proof.  The  system  (A3),  linearized  about  A,  can  be  written  as 


where  H  has  two  positive  eigenvalues  X^  ■  ^a12  “  r1  and  ^2  *  ^r2  ani1  two  n®9ative 


ones.  The  unstable  manifold  for  (A4)  Is  the  plane 


P  “  *,u,  q 


21  -  - 
171; u  -  V  • 


Its  trajectories,  projected  on  the  (u,v)  plane,  are  the  solutions  of 

-  ■  f  ■  Y  +  u  -  with  Y  “  and  w  -  IT 

p  du  u  112  1 

and  can  therefore  easily  be  investigated.  For  e  >  0,  they  cross  tranversally  the  line 


se^ent 


i  s  u”E(1  -s),  v*  -aes  (0  <  s  <  1)  , 

£  mm 


provided  that  a  >  Y.  Therefore  s  is  a  parameter  for  those  trajectories  entering  the 
quarter-plane  {u  >  0,  v  <  0}  which  cross  the  rest  cross  the  u-axis  between  0  and 

e,  and  can  be  parametrized  by  the  point  at  which  they  do  so. 

Because  of  this  transversal lty,  if  we  define  the  similar  segment  in  the  (u,v)  planet 
*e  s  u  »  e(1  -  s),  v  -  1  «  -aes  (0  <  s  <  1 ) 

for  small  enough  e,  then  s  can  be  used  to  parameterize  the  trajectories  on  M  whose 
projections  cross  t from  left  to  right  or  top  to  bottom.  We  denote  these  trajectories 
by  Xg(x)  “  (u8(x>,  ps<x),  vB<x),  q8(x)>,  and  their  projections  by  U8(x).  The  latter  are 
pictured  in  a  neighborhood  of  A  as  follows: 


(0,1-ae) 


L. 


Case  U  S  1 


(0,1-ae) 


He  are  concerned  with  elucidating  conditions  under  which,  for  some  s,  Xg  is  a 
connection  from  A  to  B.  We  shall  denote  by  Q  the  square 

Q  =  tO,  1 1  *  [0,1] 

in  the  (u,v)  plane. 

We  shall  need  to  account  for  sign  changes  of  p  and  q  along  Xg: 

Definition.  The  point  U  is  a  turning  point  for  Xg  if  it  lies  on  the  trajectory  Ug, 

and  if  ps  or  qa  vanishes  there.  (If  both  pg  and  q^  vanish  at  that  point,  it  is 

counted  twice  in  the  following  proposition.) 

Proposition  2.  If  Ug  has  two  or  more  turning  points  before  leaving  Q,  then  the  first 
two  are  well  defined  and  isolated  from  the  rest.  Immediately  after  the  second,  P8  <  0 
and  qa  >  0. 

Proof.  If  p8(x1)  -  0  for  some  Xj,  it  cannot  be  true  that  pg(x1)  *  0,  for  otherwise 
u(x)  =  u  (x„),  p(x)  =  0  would  be  a  solution  of  the  first  two  equations  in  (A3),  and  by 
uniqueness  of  the  initial  value  problem  for  that  subsystem  (v(x)  being  considered  known), 
would  equal  the  first  two  components  of  Xg( x) .  This  is  impossible  because  our 
construction  does  not  allow  u  to  be  constant.  Therefore,  p  or  q  (or  both)  changes 
sign  at  each  turning  point. 

Each  fixed  trajectory  Xg  has  p  >  0,  q  <  0  near  to  (0,1).  Suppose  p  changes 
sign  at  the  first  turning  point,  and  q  does  not.  Since  p  “  -uf  <  0  there,  that  point 
must  lie  where  f  >  0,  i.e.  below  the  line  f  ”  0.  Since  q  <  0,  the  trajectory  stays 
below  at  least  until  the  second  turning  point.  Therefore  p  cannot  change  sign  at  the 
second  point,  since  p  would  have  the  wrong  sign,  and  it  must  be  q  that  changes.  The 
other  possible  cases  are  analyzed  the  same  way.  In  particular,  although  the  first  two 
turning  points  may  coincide,  the  first  three  may  not.  This  completes  the  proof. 

Proposition  3.  If  Ug  does  not  exit  Q  and  does  not  have  more  than  one  turning  point, 
then  it  has  none,  and  Xg  is  a  connection  to  B. 

Proof.  There  being  at  most  one  turning  point,  ug{x)  and  vg(x)  are  eventually  monotone, 
so  Xg(x>  must  approach  a  rest  point  as  x  ♦  “.  There  are  four  rest  points  A,  B,  0  » 

A  A  A  A 

(0,0, 0,0)  and  E  «  (u,0,v,0)  where  (u,v)  is  the  intersection  point  of  f  «  g  »  0.  The 
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proof  consists  in  eliminating  all  except  B.  A  is  excluded  because  it  would  require  two 
turning  points.  We  show  why  E  is  also  excluded:  the  argument  for  0  is  similar.  If 
t!a  ♦  (u,v)  as  x  ♦  the  approach  must  be  monotone.  It  cannot  be  from  the  left  and 
above  f  ■  0,  because  then  p  >  0  and 

p  «  -uf  >  0,  contradicting  the  requirement  that  p  ♦  0.  Similarly  it  can  not  approach 
from  above,  but  remaining  below  g  •  0.  This  excludes  all  approaches  from  above  left. 
Approaches  from  below  left  are  impossible  because  on  the  approach,  we  would  have  q  >  0. 
Since  only  one  sign  change  was  allowed,  it  must  have  been  q,  and  the  turning  point  must 
have  been  below  f  =  0  and  g  »  0.  But  that  is  impossible  because  q  would  have  the 
wrong  sign  there.  The  other  directions  of  approach  may  be  eliminated  for  similar 
reasons.  This  completes  the  proof. 

Definition:  Oq  e  a  is  the  class  of  functions  (f,g)  such  that  for  some  s  6  (0,1),  Ug 
has  two  or  more  turning  points  before  leaving  Q. 

Proposition  4.  Suppose  (f,g)  ^  0'  and  there  no  monotone  connection,  A  to  B.  Then 

every  (J8  exits  Q  transver sally . 

Proof.  By  Proposition  3,  Ug  must  exit  0  before  its  second  turning  point  (if  it  has 
one).  If  it  exits  through  the  Interior  of  the  bottom  of  0,  then  q  <  0  just  before  and 
after  exit*  if  q  -  0  at  exit,  then  q  -  0  there,  and  the  last  two  equations  of  (A3) 
would  attain  a  rest  point  at  exit,  which  is  impossible.  Therefore  q  <  0.  Similarly  exit 
through  the  interior  of  the  right  side  can  not  be  achieved  unless  p  >  0  there.  The  other 
cases  are  handled  the  same  way. 

Proposition  5.  a  0  is  an  open  set  in  the  space  of  parameters  rA  and  a^  (i,j  “  1,2). 
Proof :  Consider  a  trajectory  Xg  with  two  turning  points.  The  proof  of  Proposition  4 
also  shows  that  sign  change  of  pg  can  not  occur  on  the  lateral  sides  of  Q,  nor  those 
of  qg  on  the  top  or  bottom.  If  pg  changes  sign  before  qg,  it  must  therefore  do  so  in 
the  Interior  of  Q,  and  then  the  second  turning  point  must  be  in  the  interior  or  on  the 
left  side.  The  latter  is  excluded  because  qg  -  -v^g  <  0  there.  A  similar  argument  holds 
if  the  order  is  reversed.  Therefore,  both  turning  points  occur  in  the  interior.  Mow  it  is 


easy  to  see  that  the  trajectory  Xg,  with  s  fixed,  depends  continuously  on  the 
parameters  of  the  equation.  It  therefore  continues  ho  have  two  turning  points  before 
exiting  p  when  the  parameters  are  changed  by  a  small  amount. 

Theorem  A 1 .  Let  the  family  (f,g>r  e  d  depend  continuously  on  the  real  parameter  r. 

Assume  (f«g)r  *  Oq  for  r  »  0,  but  not  for  all  positive  r.  Let  r*  = 

sup{ r  :  (f,g)  e  a„  for  0  <  r  <  r} .  Then  A  is  connected  to  B  for  r  =  r*. 

x  0  *  * 

Proof.  Our  trajectories  now  depend  on  two  parameters,  that  is,  Ug  r*  Let  {r^}  *>e  a 

sequence  with  r  t  r*,  and  let  s  be  such  that  U  has  two  turning  points  before 

"  n  sn'rn 

leaving  Q.  By  possibly  taking  a  subsequence,  we  may  assume  that  s^  ♦  s*  as  n  ♦  ». 

Now  a*  e  (0,1),  because  it  is  easily  seen  that  any  trajectory  U  with  s  in  a  neigh- 

borhood  of  0  or  of  1  does  not  have  two  turning  points  before  leaving  Q.  Consider  the 

trajectory  Xa#  r#.  By  Proposition  5,  (f,g)r*  ^  hence  by  Proposition  4,  either  there 

exists  a  connection  for  r  =  r* ,  or  every  trajectory  U  *  leaves  0  transver sally , 

before  experiencing  two  turning  points.  The  proof  of  the  theorem  consists  in  excluding  the 

second  case.  That  case  would  imply,  in  particular,  that  0  *  exits  transversally ,  and 

that  while  us*#r*(x)  C  Q,  either  p  or  q  remains  of  one  sign  and  is  bounded  away  from 

zero  near  the  exit  point.  By  continuity  of  trajectories  with  respect  to  parameters, 

X  .  #  can  be  approximated  by  portions  of  the  trajectories  X  for  large  n. 

■  ,r  n'n 

Therefore  for  such  n,  the  also  exit  transversally  before  experiencing  a  second 

turning  point.  This  contradicts  their  construction  and  finishes  the  proof. 

Corollary  to  Proposition  4 .  If  (f,g)  ^  0^,  then  there  exists  a  monotone  trajectory  which 
attains  or  approaches  (1,0).  If,  in  addition,  no  connection  exists,  then  any  given  point 
on  the  bottom  or  right  side  of  Q  is  attained  by  some  trajectory  Ug. 

Proof ;  If  there  exists  a  monotone  connection,  we  are  done.  If  not.  Proposition  4  shows 
that  every  Ua  exits  transversally.  For  small  s,  the  exit  point  is  on  the  top  of  Q. 

For  s  “  1  the  trajectory  lies  on  the  v  axis  and  the  exit  point  is  the  origin.  By 
transversality,  a  standard  argument  shows  that  the  exit  point  depends  continuously  on  s. 
Therefore  the  range  of  exit  points  must  include  every  point  between  those  for  s  near  0 


and  s  near  1,  proceeding  either  clockwise  or  counter.  But  (0,1)  is  not  an  exit  point 
(it  would  require  two  turning  points),  so  we  proceed  clockwise,  and  cover  all  of  the  sides 


indicated. 

For  the  following  proposition,  we  divide  Q  into  four  regions 

I  -  {f  <  0  <  g}  n  0,  II  *  (f  >  0,  g  >  0}  0  Q  , 

III  -  {g  <  0  <  f)  fl  Q,  IV  »  (f  <  0,  g  <  0}  n  8  . 

For  a  monotone  trajectory  Ug  from  (0,1)  to  (1,0),  we  denote 

r  =  {u  <x>}  n  c, 

s 

JIf  “  I / m j “f < u » v ) du | ,  J  -  l/rnIvg(u,v)dv| 


and  similarly  for  Jixf'  etc* 

Proposition  6.  Let  there  exist  a  monotone  trajectory  Ug  approaching  or  attaining 
(1,0).  Then  the  following  estimates  hold,  where  the  symbols  c  and  C  denote  positive 
constants  depending  only  on  lower  and  upper  bounds  for  the  constants  r^,  a^  (i,j  ”  1,2), 
and  a,  where  we  define  a  ■  the  length  of  dill  (v”0 } .  They  do  not  represent  the  same 
constants  in  each  Instance.  Here  b  is  the  diameter  of  the  region  I. 


(a)  Jlf  ♦  JIg  <  Cb> 


(b)  Jjyg  *  c^jvf  ™ 


(c)  the  inequality  +  JIVf  -  8  <  °  implies 


JIIIg  +  JIVg  >  c  and  JIIg  i  c(s+b)' 


(d)  JJf  +  Jivf  >  JJIf  +  J1IIf) 


,e>  Jig  +  Jug  >  Juig  +  Jivg- 


Proof ;  (a)  follows  immediately  from  the  fact  that  the  length  of  T  fl  I  is  bounded  by  2b. 

(b)  In  the  limiting  case  I  “  9,  we  have  vg  >  cuf  in  IV,  from  which  the  result 
follows  in  that  case  with  C  “  0.  When  1/9,  this  inequality  is  violated  only  near  the 
intersection  91  fl  9111,  which  has  length  <  b. 

(c)  The  following  sketched  argument  can  easily  be  quantified.  The  integral 


•^ro(rituiV) ,uf  Idu  “  Jnif  *  Jivf  can  1,6  made  8Bla11  only  by  requiring  that  T  O  (III  V)  IV) 


lie  mostly  in  some  small  enough  neighborhood  of  {f  *  0),  the  sire  of  the  neighborhood 


depending  on  a  lower  bound  for  a.  But  lying  in  this  neighborhood  forces  I /  vgdv|  to  be 
larger  than  some  minimum  amount,  since  |vg|  is  not  also  small  in  that  neighborhood.  It, 
together  with  the  monotoncity  of  V,  also  forces  T  ft  II  to  lie  near  9 1  (ft  9 II,  hence 
the  final  estimate  in  (c). 

(d)  Multiply  the  second  equation  of  (A3)  by  p  and  integrate  from  -•  until  (1,0) 
is  attained.  We  find  that  the  value  of  p  at  (1,0)  satisfies 

p2  =  ~/r  uf (u,v)du  >  0  . 

The  integral  appearing  here  is  equal  to 

"  JIf  “  JIVf  +  JIIf  +  JIIIf  ' 

hence  (d). 

(e)  is  proved  in  the  same  manner. 

For  purpose  of  the  following  theorem,  we  remove  the  restriction  *2^*22  “  r/*n  “  *• 
Proposition  6  is  clearly  valid  without  it. 

Theorem  A2 .  Let  the  coefficients  rt  >  0  and  a^  >  0  (i,j  “  1,2)  depend  continuously  on 
a  real  parameter  £  for  £  in  some  interval  A.  Assume  that  for  all  £  e  A 

a22  .  Vi 

a12  *11 

Assume  that  for  some  £  e  A,  <  *22^*12'  an<J  f°r  another  value»  r2^ri  »  a21/all’ 

Then  for  some  £  e  A,  there  exists  a  solution  of  (Al),  (A2). 

Proof:  For  values  of  £  such  that  r 2/r  >  *21^*11'  we  have  that  (f,g)  e  a.  but  ^  «q. 
To  see  the  latter,  observe  from  the  proof  of  Proposition  2  that  if  a  trajectory  has  two 
turning  points,  they  must  lie  in  the  region  III  where  f  >  0  and  g  <  0.  But  the 
assumed  inequalities  imply  that  region  is  empty. 

On  the  other  hand,  by  continuous  dependence  on  £,  there  are  values  of  £  for 
which  a22^a12  <  r2/,r1  <  a2 i^ai  1  and  r2/,r1  is  arbitrarily  close  to  a22/a21.  When  it  is 
close,  the  size  of  I  is  small,  the  parameter  a,  remaining  bounded  away  from  zero.  For 
sufficiently  small  b,  we  show  (f,g)  e  aQ.  This,  together  with  Theorem  Al  will  establish 
the  desired  result.  It  suffices  to  show  there  exists  no  monotone  trajectory  attaining  or 
approaching  (r^/a.,^,0).  Suppose  there  were.  A  straight-forward  calculation  shows  that 


r  must  be  concave  upward  in  IV  and  downward  in  II.  It  follows  from  this  and  the  fact  that 


T  must  reach  (r^/a^,*))  that  either  T  II  or  T  fi  IV  must  be  empty.  Suppose  it  is 
the  former.  For  some  S  >  0,  let  be  the  set  of  C  for  which 

*22^*12  <  r2^r 1  <  *22^*12  +  In  the  following,  the  constants  c  and  C  are  uniformly 

bounded  from  below  and  above  for  C  e  A  ,  and  do  not  necessarily  denote  the  same  constant, 

o 

even  within  a  single  context.  From  (a)  in  Proposition  6,  we  have  Jjj  ♦  JIg  <  Cb,  and 
from  (e). 


Cb  >  J  >  JTrT  +  J 

Ig  Illg  IVg 


Hence  if  b  is  small  enough,  the  inequality  in  (c)  regarding  +  Jxvg  is  violated. 

If  +  JIVf  *  Juif#  we  have 


hence  JIJIf  ♦  JjVf  *  c*  Fron  this  and  the  fact  (d)  that  JTJ!  +  J 


j  +2J  >  j  J  >c 

If  IVf  Illf  IVf 


JIVf  >  2  (C  ‘  <*>  ‘ 

c  2 

The  same  type  of  inequality  also  holds  (b)  for  JIVg  s  JIVg  *  J  (c  “  Cb)  >  c  /4  (for  b 

2 

small  enough).  From  (e),  this  implies  J  >  c  /4,  which  for  small  b  contradicts  the 

ig 

fact  that  J_  <  Cb. 

ig 

The  other  case  is  when  r  fl  IV  -  $.  Again,  we  have  <  Cb.  From  (d),  we  obtain 

Cb  >  Jif  >  JIIf  "  Jjxf  +  Jxvf '  80  the  hypothesis  of  (c)  is  fulfilled  for  small  b.  Hence 

J  <  C(b  +  s)  <  Cb  . 

Hg 

But  this  contradicts  (e),  whose  left  side  is  bounded  by  Cb,  and  right  side  (by  (c) )  is 
bounded  from  below.  Therefore,  no  such  T  can  exist.  This  completes  the  proof. 
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